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Abstract 

In [6], we proved an asymptotic 0(n -Q// ( a+1 )) bound for the approximation 
of SU(./V) loops (N > 2) with Lipschitz smoothness a > 1/2 by polynomial 
loops of degree < n. The proof combined factorizations of SU(iV) loops into 
products of constant SU(iV) matrices and loops of the form e A ^ where A{t) 
are essentially su(2) loops preserving the Lipschitz smoothness, and the careful 
estimation of errors induced by approximating matrix exponentials by first- 
order splitting methods. In the present note we show that using higher order 
splitting methods allows us to improve the initial estimates from |6| to close- 
to-optimal ©(ra - ^ -6 )) bounds for a > 1, where e > can be chosen arbitrarily 
small. 



1 Introduction 

The study of approximation rates for Lie-group-valued loops by polynomial loops 
is a relatively unexplored topic within the larger area of nonlinearly constrained 
approximation. Motivation is provided by previous density results [21 El E] for semi- 
simple Lie groups, and by more practical needs, e.g., for the design of para-unitary 
FIR filters [21 El 10] . 

In this note, we continue the study of the SU(AQ case, N > 2, and improve 
upon the following Jackson-type estimate stated in |BJ: : For any Lip a - continuous 
loop U : TT — > SU(iV) and a > 1/2 there exists a sequence of polynomial loops 
U n : TT — > SU(A^) of degree < n such that the following asymptotic inequality holds 

\\U - U n \\ c := max \\U(t) - U n (t)\\ < C ajNtU (n + n >0. (1) 

Even though this estimate is admittedly far from final, to our knowledge it repre- 
sented the first nontrivial upper estimate for the achievable rate of approximation 
for Lip a loops with values in matrix Lie groups. Note that there is a large gap 
between the exponent a /(a + 1) established in (OQ), and the trivial upper bound a 
for the maximal order of approximation of Lip a loops following from the classical 
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Jackson-Bernstein theorems for univariate trigonometric approximation [TJ Theorem 
7.3.3]. 

In the meantime, we realized that some simple modifications in the proof strategy 
of [6J yield close-to-optimal rates, at least for a > 1. The major change is to use 
higher order splitting methods for exponentials instead of the standard first order 
approximation e^ =1 ^ « e Xl . . . e Xj . 

Theorem 1 Let N > 2, a > 1, and e > 0. For any U(t) G LipjT -> SU(iV)), 
there exists a sequence of SU(iV)- valued polynomial loops U n (t) of degree < n such 
that 

\\U - U n \\ c < C a;£iN;U (n + n > 0. (2) 

The proof of this result is given below. The restriction a > 1 comes from the 
fact that we currently miss error formulas for splitting methods of order k > 2 in 
terms of higher order commutators. Another obstacle is the lack of formal proof for 
factorizations of SU(iV)-valued Lip a loops into exponentials of su(iV)-valued Lip a 
loops if a < 1/2. The latter problem can possibly be removed by using homotopy 
arguments as in [3] (the first author acknowledges inspiring discussions with W. M. 
Lawton on this and related subjects of the present note). See also the remarks at 
the end of the next section. 

The major open question is whether (j2J) remains true also with e = 0, or if 
the nonlinear constraints lead to a slight deterioration of approximation results 
compared to the unconstrained case. Settling this question will probably require a 
different approach. 



2 Proof of Theorem 1 

We first recall the facts already proved in [6]. The notation we use is either self- 
explaining or can be found in [S] (we have opted to keep notation very close to that 
of [B], to make the comparison easy). The matrix norm of choice is the spectral 
norm. 

Factorization into essentially exponentials of su(2) loops. Lemma 4 of [6J 
states that for any a > 1/2 and any U(t) e Lip Q (Tr — ► SU(iV)), there exist constant 
matrices U 0> i G SU(iV) and loops A t (t) G Lip a (T -> su(2)) such that 

L 

U(t) = n Uwe* 1 ®, teT, L := N(N - l)/2, (3) 
i=i 

where A\(t) = TijAi(t) denotes the canonical extension of Ai(t) to a su(iV) loop by 
the map 
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for some index pair with 1 < % < j < N (If. denotes the k x k identity matrix). 
Moreover, smoothness of the factors is controlled by smoothness of U(t): 

||A|kip tt <C ai7V)[ ;||[/||Lip a) l = l,...,L. (4) 

Approximation can be done factor-by- factor. If Pi(t) are polynomial loops 
in SU(2) of degree < n such that 

\\e Mt) - < l = l,...,L, (5) 

for the Ai(t) occurring in the factorization ([3]) then 

L 

p(t) = n u oAv)> ( Pi® = T ij p i(t) ) 
i=i 

is a polynomial loop in SU(iV) of degree < Ln and it satisfies the estimate 

\\U(t) - P(t)\\ c < Le. (6) 

Use Lemma 5 from [B]. 

Construction of Pi(t). For any m > 1, we can approximate Ai(t) G Lip^ClT — > 
su(2)) by a su(2)-valued polynomial loop Ri, m (t) of degree < m at optimal rate (say, 
by applying Vallee-Poussin means componentwise), i.e., 

|| e 4(0 _ e^-^Hc < Cm"", m > 0, (7) 

see Lemma 6 a) in [6j, and 

6 m 

Rl,m(t) = Cr,kBr,k(t) 
r=l k=0 

(to keep notation simple, the dependence of the coefficients Cfc jr on I and m is not 
made explicit, moreover, the terms with k = and r even are fake). Here {S rj fc(t)} 
is the designated basis (over R) for the linear space of dimension 6m + 3 of all su(2)- 
valued polynomial loops of degree < m, see [6] . Lemmas 1 and 4 in [6] establish the 
following facts which are relevant below: 

\\Cr, k B r>k (t)\\c = M <C(fc + l)- Q , k = 0,...,m, r = l,...,6, (8) 

and if ordered properly (e.g., lexicographically), products of the form 

6 m 

c rk B r ^(t) 
r=l k=0 

represent SU(2)-valued polynomial loops of degree < 6m, independently of the choice 
of A > and the set of real coefficients {c r> k}- For our applications we take A = 1/M, 
for some integer M > 1. 
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With these preparations, we can write down the final formula for the approxi- 
mation of e Al ^: 

Pi(t) :=0({c r , fcJ B rifc (t)/M} r=li ... i6;fc=o ,..., m ) M , (9) 

where <f>{{Xj}j=i,...,j) is a suitable splitting method for approximating e Xi+—+Xj^ 
see the next paragraph for details. The integers m and M will be fixed later. Note 
that the pointwise estimate 

\\ e Mt) _ m \\ < Cm ~ a + M\\e R ^ M - 0({c r , fcJ B r , fc (t)/M} r=1 ,..., 6;fc=o ,..., m )|| (10) 

follows from applying the triangle inequality to 

e Mt) _ p^) = ( e Mt) _ e R l<m (t)j + ^ e R l>m (t)/MjM _ 

then using (|7j) for the first term, and Lemma 5 in [6] for the second. Thus, the 
quality of approximation crucially depends on the properties of the chosen splitting 
method <fi. 

Estimate for the second term in (1 10ft . Now we depart from [6J, where the 
method <p of choice was the first order splitting method 

o.i{.Y,i, , /) : r v V v -'... ( v <. 

The error estimate is stated in Lemma 6 b) of [U], and leads to an overall estimate 
< CM -1 for the second term in the right-hand side of ( jTUl) if a > 1/2 (and to the 
suboptimal asymptotic approximation rate of that paper). We now show that using 
higher order symmetric methods leads to significant improvements. The standard 
2nd order symmetric method is given by 



Xl/2 . . . e Xj - l/2 e Xj e Xj ~ l/2 . . . e Xl/2 



M{ x j}k=i,...,j) ■= e 

Following Yoshida (see [1]), symmetric splitting methods of order 2(s + 1) can be 
constructed from a given symmetric method of order 2s via the formula 

02(s+i)({^j}j=i,...,j) := 02s({a^}i=i,...,j)02s({&sX i } i= i i ... i j)02,({a s X i } i= i i ... i j), 

if one chooses the constants as follows: 

a s = (2 - 2 1 /( 2 *+ 1 ))-i, b s = -2 1 /( 2s +D(2 - 2 1 /^)- 1 . 

The order condition for these (f>2s can be stated as follows: For A — > 0, we have 

l|e A(Xl+ - + ^ ) -^({X J } J=1 ,..., J )|| = 0(A 2 ^). 

Using Taylor expansion and rough estimates, the order requirement translates into 
the error bound 

2s+l 

A(x 1+ ...+x J ) _ fafixy < CA 2S+1 [ V IIX-II 1 , (ID 



^({x,} j=w )||<ca 2s+i (em) 



valid for |A| < A max , with a constant depending on A max and J?. \\Xj\\. More precise 
error bounds are available for 2 > see Remark 2 at the end. 

We are now ready to apply this to the family {c rj jfcS r . jfe (t)} r=li _ )6 . k=o,...,m (point- 
wise with respect to t) with A = M" 1 . By ([H]), if a > 1 we have 

6 m 

Y^Y,W C ^ B rAt)\\<C, (12) 
r=l fc=0 

where C depends on the Lip a -norms of the the su(2)-loops Ai(t), but not on m. 
Substituting into ( TTTT) . we obtain for each I = 1, . . . , L 

\\ e R, m (t)/M _ 2a ({ Cr)fc S r)fc (i)/M} r= i i ... )6 . k=o,...,n) lie < CM-< 2s+1 >. (13) 

From now on, set cf) = fas in the formula for Pi(t) (and consequently P(t)). Substi- 
tuting into (fTUl) and taking into account © we finally arrive at 

\\U{t) - P{t)\\ c < Lmax \\e Mt) - P z (*)|| < C(m- Q + M" 2s ), (14) 

where C depends on s, a, iV, and on U(t). 

Estimating the degree of P(t). Consider a large enough integer n > hq 
(for n < no, just use constant P(t) = I to get the complementing trivial bound 
\\U(t) — I\\ < 2). We will now fix m and M such that the degree of the above 
constructed P(t) is < n and the right-hand side in (jl4jl is asymptotically as small 
as possible. Due to the recursive definition of fas, the degree of the polynomial 
loops (p2s({c r ,kB r ,k{t) I M}) is bounded by 3 s " 1 times the degree of the polynomial 
loops (p2{{c r ,kBr,k{t) / M}) generated by the 2nd order method (for simplicity, we do 
not indicate the index set r = 1, . . . , 6; k — 0, . . . , m in the notation). The latter, 
however, have degree < 12m. This can be proved as in Lemma 1 of [6J. Indeed, 
write 

fo({cr,kB r , k (t)/M}) = 1 ({c P>fc S rjfe (t)/2M})0 1 ({c r , fc S Pifc (t)72Af})*. 

We already know that the first factor <fii({c r ^B r ^(t) /2M}) has degree < 6m. The 
second factor is the Hermitean transpose of (f>i{{c r ^B* k (t) /2M}) , and it remains to 
check that {B* k (t)} is such a permutation of the original basis {B rt k(t)} to which 
Lemma 1 of [6] can be applied, leading to the same degree bound. 

Putting things together, we see that the degree of P(t) is bounded by 12L3 s_1 Mm. 
Thus, choosing the integers M, m according to 

M = [(12L3'- 1 )- 1 n a/(a+28) ], m = [n 2s ^ a+2 % 

we guarantee that the degree of Pit) does not exceed n. On the other hand, substi- 
tution into (1141) yields 

\\U(t) - P(t)\\ c < Cn~ 2sa ^ a+2 ^ = Cn- a+a2 ^ a+2s \ 

This establishes the claim of our theorem, if, for given a > 1 and e > 0, we choose 
the order 2s of the splitting method large enough. 
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Remark 1. There are at least three shortcomings of the asymptotic estimate 
(j2J). First, the constant C a>€j N,u depends on U(t) in an unspecified way. Secondly 
the restriction a > 1 is mainly due to the use of the crude error estimate fill I) for 
higher order splitting methods (see the comments in Remark 2). In addition, for 
a < 1/2 the factorization technique of Lemma 4 from [6 J breaks down (an alternative 
is addressed in Remark 3). Thirdly, it is not clear at the moment if one can set e = 
in ©• 

Remark 2. For low-order splitting methods such as 0i and 02, the error bound 
can be made more precise in terms of commutator expressions which paves the way 
for proving (1131) (and thus also (|T4"|) ) for some a < 1. For 0! this was demonstrated 
in [Hj (see Lemma 6 for the more precise error bounds). An improved error bound 
for the symmetric 2nd order method 02 has been established in [9]: 

|| e A(x 1+ ... + x J ) _0 2 ({^.} i=w )|| < A 3 A(Xi, . . . ,Xj), (15) 
where A{X 1 , ...,Xj)= Yf k Z\ A 2 (X fc , X k+l + ■ • ■ + Xj), and 

A 2 (A,B) = ±{\\[[A,B],B]\\ + ±\\[[A,B],A]\\}. 

The advantage is that in our application, where we would apply ffToT) to matrix sets 
of the form {ckB r k}k =Qt . m related to the terms of the Fourier series of a su(2)- 
valued Lip a -loop, the sum of the norms of the appearing triple commutators can be 
estimated by a sum of the form X^fcLiO°g k) 2 k~ 3a which remains uniformly bounded 
for a > 1/3 (the details are worked out in [6] for the first-order case). In contrast, 
using fTTTj) with s = 1 leads to a constant factor of the form (J2T=i k~ a ) 3 - 

Unfortunately, generalizations of (j!5p to Yoshida-type (and any other higher- 
order) splitting methods are not known. 

Remark 3. W. R. Lawton drew our attention to a possible alternative to loop 
factorizations of the form (J3J) proposed in [S]. It is well known that SU(iV) is a 
simply connected compact C°°-manifold. Thus, any SU(iV)-valued continuous loop 
U(t) can be contracted to a point by a homotopy map ip : [0, 1] — > C(TT — > SU(iV)) 
(i.e., ip is continuous, -0(1) — U(t), and (without loss of generality) -0(0) = /)• 
Let us assume that for U{t) G Lip^TT — > SU(iV)), the homotopy map ip can be 
found in such a way that ■?/>(£)€ Lip a (TT — > SU(iV)) for all £ G [0, 1] (i.e., preserves 
Lipschitz smoothness along the homotopy path). We do not have a reference for 
this assumption but strongly believe that it holds for all a > 0. 

Now, take a fine enough partition £ = < £i < ... < £k-i < £k — 1 of 
[0, 1] such that H^^-i) — ip(£,k)\\c < r N, where r N is the injectivity radius of the 
exponential map in the neighborhood of / G SU(iV). Then we can write 

U(t) = U t (t) . . . U K (t), U k (t) := *l>(tk-i)*il>(Zk), k = 1, . . . , K, 
where all Uk(t) belong to Lip^H 1 — > SU(iV)), and 

||/ - U k \\ c < HVfe-irilcllV'fe-i - ^k)\\c = - 4>(£k)\\c < r N , 
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i.e., A k (t) = log([/fc(t)) is well defined and belongs to Lip a ( n D n — > su(iV)) for all 
k = 1, . . . , K. Thus, 

K 

U{t) = \[e A ^\ \\A k \\ LiPa <C(a,N,U), k = l,...,K. (16) 

k=l 

In contrast to ([3]), the number of exponential factors K is not independent of a and 
U, and the A k (t) are general su(iV)-valued Lip Q -loops (and not essentially su(2)- 
valued as in Q). 

However, as long as we accept the dependence on U{t) in the constant appearing 
in (121), the factorization ( fTBT) is sufficient to carry out the above proof with minor 
changes. The reduction to the SU(2) case can be circumvented by working with a 
similar basis {-Br,fc(^)}r=i,...,_Rjv,fc=o,...,m over R for SU(iV)-valued polynomial loops of 
degree < m. What changes is the number R N of subsets {-B r ,fc(^)}fc=o,...,m of basis 
elements to be considered. This number depends only on N, and enters the degree 
estimates as a linear factor. 

Remark 4. Loop approximation can be pursued for other Lie groups and man- 
ifolds as well. Work on the closely related case of SO(iV)-valued loops (N > 3) is 
ongoing [H]. 
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